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Introduction: 

DSMC  Chemistry  Models 


•  Important  for  hypersonic  reacting  flows  (re-entry):  strong 
nonequilibrium,  large  effect  on  heating,  radiation 

•  Long  history  of  proposed  simplified  models;  realism  can 
be  improved  if  reaction  behavior  is  well  understood 

•  After  flurry  of  work  inspired  by  BSUV  data,  many  felt  that 
without  new  measurements,  further  work  unfruitful 

•  New(ish)  Quantum  Kinetic  (QK)  model  proposed  by  Bird 
has  generated  renewed  interest 
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QK  Model 


*  See  G.A.  Bird,  Gallis,  references  therein 

*  Very  appealing  simplicity 

*  Impressively  predicts  thermal  rate  coefficients  k(T) 
over  wide  T  range  in  good  agreement  with 
measured/recommended  values 

*  The  present  paper  investigates  the  behavior  of  o(E), 
or  reaction  probability  per  collision,  important  for 
strongly  nonequilibrium  flows 

*  Available  detailed  quasi-classical  trajectory  (QCT) 
results  for  o(E)  for  2  important  reactions  are  used  to 
evaluate 
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Model  Definitions 


*  QK  dissociation:  Preact  =1  for  Ec>Ea,  =0  otherwise 

*  QK  exchange: 

*  Total  Collision  Energy  (TCE)*:  reproduces 

—  k(T)  =  ATBexp(-Ea/T)  for  given  values  of  A,  B,  Ea  and  C, 

*  Both  models  need  total  collision  cross  section 
values:  VHS  model  parameters  selected  so  that 
collision-pair  viscosity  matches  values  predicted  by 
Capitelli  (ref)  in  5,000  -  15,000K  range 

*  *  Bird,  1994. 
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Vibrational  Energy  Levels 


*  Original  TCE  model  used  continuous  internal  energy 
modes 

*  Later  work  adapted  it  for  quantized  internal  energy 
levels,  see  esp  Gimelshein  (ref)  and  discussion  of 
small  correction  factor  to  use  for  simple  harmonic 
oscillator  (SHO)  levels;  zero-point  vibrational  energy 
(ZPE)  is  included  in  total  collision  energy 

*  QCT  calculations  of  course  use  real  vibrational 
energy  level  values  (closer  to  AHO) 

*  This  work  compares  DSMC  results  for  SHO  levels  to 
QCT  results  for  AHO  levels  which  have  same  Evib 
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N2+N  Dissociation  Reaction 


*  Plot  of  thermal  reaction  rate  coefficients 

*  TCE  model  uses  published  Arrhenius  parameters  as 
input  and  will  match  QCT  k(T)  by  definition 

*  QK  model  uses  no  adjustable  parameters  and 
matches  QCT  (and  expt)  k(T)  quite  closely 


Cross  section,  m 


N2+N  Dissociation  Reaction 


*  QCT  o  results  from  Esposito  (Chem  Phys  Lett  302, 1999  and  unpublished) 

*  Dissociation  reaction,  as  is  well-known,  has  strong 


vibrational  favoring 


8 


N2+N  Dissociation  Reaction 


*  State-specific  thermal  rates  kv(T) 

*  TCE  and  QK  models  both  have  kv(T)  that  are  high  for 
low-v  and  low  for  high-v  levels.  No  surprise. 
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N2+N  Dissociation  Reaction 


QK  assumes  that  each 
collision  above  threshold 
energy  has  unity  reaction 
probability 

TCE  matches  realistic  QCT 
o  surprisingly  closely 
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N2+N  Dissociation  Reaction 


For  v=60,  QCT  has  odiss  ^ 
ototai ;  QK  very  close  to  this 

TCE  far  underestimates  the 
probability  for  very  high  v  . 


For  flows  with  mainly 
low  v  populated,  QK 
may  be  problematic 
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N  +0  Exchange  Reaction  k(T) 


*  N2+0  ->  NO  +  N  QCT  results  from  Bose  and  Candler, 
unpublished 

*  Plot  of  thermal  reaction  rate  coefficients 

*  TCE  model  uses  published  Arrhenius  parameters  as  input  and 
will  match  QCT  k(T)  by  definition 

*  QK  model  uses  no  adjustable  parameters  and  matches  QCT 

(and  expt)  k(T)  reasonably  — , — — , — — , — — , — — , — 


-  Bose  refit 
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Cross  section,  cm 


N2+0  Exchange  Reaction 


*  This  (and  many  other)  exchange  reactions  have  no 
strong  vibrational  favoring 
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QCT  cross  sections  at 
same  Ec  show  ~3x 
difference  in  o 

Lends  itself  more 
readily  to  a  simplified 
DSMC  model 


Approved  for  public  release,  distribution  is  unlimited 


13 


N2+0  Exchange  Reaction 


*  Overall  shape  of  o(E)  much  better  represented  by 
TCE  model 
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Cross  Section, 


N2+0  Exchange  Reaction 


*  Overall  shape  of  o(E)  much  better  represented  by 
TCE  model,  v=0  and  v=5  examples  (J=0) 
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Sensitivity  of  Flowfield  Results  to 
N2+0  reaction  model 


*  Simple  axisymmetric  flow  over  cylinder,  8  km/s,  120 
km  (fractions  of  N2,  O?)  nearly  single-collision 


*  QCT  cross  sections  were  used  in  one  run  as  look-up 
table  with  interpolation  for  comparison 


Sensitivity  of  Flowfield  Results 


*  NO  product  mole  fraction,  8  km/s 
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•Conclusion:  As  indicated 
by  o  plot,  QK  model 
significantly  overpredicts 
reaction  probability  near 
threshold 

•QK  NO  fraction  is  3x  high 
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Sensitivity  of  Flowfield  Results  to 
N2+Q  reaction  model 

*  Simple  axisymmetric  flow  over  cylinder,  6  km/s,  120 
km  (fractions  of  N2,  O?)  nearly  single-collision 
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Sensitivity  of  Flowfield  Results 


*  NO  product  mole  fraction,  6  km/s 
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Conclusion:  As 
indicated  by  o  plot,  QK 
model  significantly 
overpredicts  reaction 
probability  near 
threshold 
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Summary  and  Recommendations 


*  Present  results  are  complementary  to  previous 
publications  on  QK  model,  focusing  on  noneq 
aspects 

*  for  strongly  nonequilibrium  flows,  the  probability  of 
reaction  is  sensitive  to  the  shape  of  the  energy- 
dependent  cross  section;  TCE  more  realistic  than 
QK  model  in  this  respect 

*  Possibly  explore  gaining  the  best  of  both? 
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